Defect induced rigidity enhancement in layered semiconductors 
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Abstract 

We discuss the mechanism responsible for the observed improvement in the structural properties 
of In doped GaSe, a layered material of great current interest. Formation energy calculations show 
that by tuning the Fermi energy. In can substitute for Ga or can go as an interstitial charged 
defect ^Inj*"^ j . We find that In?"*" dramatically increases the shear stiffness of GaSe, explaining the 
observed enhancement in the rigidity of In doped p-GaSe. The mechanism responsible for rigidity 
enhancement discussed here is quite general and applicable to a large class of layered solids with 
weak interlayer bonding. 

PACS numbers: 61.72.Bb, 62.20.de, 61.50.Ah, 71.55.Ht 
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Defects can dramatically alter the physical properties of solids in general, semiconductors 
in particular. Whereas the electronic transport properties of semiconductors are controlled 
by defects (doping and scattering), their mechanical properties are less so. In this letter 
we show that the mechanical properties of GaSe, a layered semiconductor of great current 
interest, can be controlled by suitable choice of defects and the Fermi energy. The idea is 
very general and can be applied to a large class of soft layered solids. 

GaSe is one of the members of a class of semiconducting compounds formed out of group 
III (Al, Ga, In) and group VI (Se, Te) elements. These are quasi-two dimensional lay- 
ered systems containing blocks consisting of four atomic planes (Se-Ga-Ga-Se). The atoms 
inside a block are strongly bonded (covalent and ionic) whereas the interaction between 
different blocks is weak Van der Waals (VdW) type. GaSe is a highly efficient nonlinear 
optical material, with applications in second harmonic generation, frequency mixing, and 
generation/detection of terahertz radiation . However due to the weak VdW interac- 
tion between the blocks its mechanical properties such as hardness and cleavability are not 
satisfactory for practical applications. It has nearly zero hardness by Mosh scale, cleaves 
easily along the planes parallel to the atomic layers and the nonlinear properties are difficult 
to reproduce from sample to sample [1]. These difficulties hamper the use of large GaSe 
crystals in practical applications. Dramatic improvement in the crystal quality has been 
reported in GaSe doped with 0.1 — 3 mass% In j^. Since In and Ga are isovalent, why In 
doping increases the rigidity of GaSe is a big puzzle. In this letter we discuss the results of 
ab initio electronic structure calculations involving In defects, explaining this puzzle. The 
fundamental question is where do the In impurities go in the host lattice and whether or not 
they strongly affect the local bonding between both host-host and host-defect pairs leading 
to the hardening of the elastic constants. We show that by a suitable control of the location 
and the charge state of the In defect one can indeed dramatically enhance the interlayer 
rigidity of GaSe. 

Due to its layered structure GaSe crystallizes in different polytypes. Since the physical 
properties of the different polytypes are quite similar we have chosen f3 — GaSe as the test 
case. /3— GaSe crystallizes in layered, hexagonal structure, having the space group PG^/mmc. 
As mentioned before GaSe consist of blocks of four atomic planes; two planes of Ga atoms 
sandwiched between two planes of Se atoms as illustrated in Fig. [T](a). In each plane the Ga 
or Se atoms are arranged in a two-dimensional hexagonal lattice. An important feature of 
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this compound is the existence of Ga-Ga dimers oriented perpendicular to the layers; their 
presence leads to different local geometry of Ga and Se. Ga has three Se nearest neighbors 
(NN) and one Ga NN, while each Se is bonded to three Ga atoms and no Se atoms. The 
presence of Ga-Ga dimers leads to bonding and antibonding Ga 4s states. The antibonding 
4s strongly hybridizes with the Se p bands and opens up a gap near the Fermi energy, leading 
to a semiconductor [4]. It is interesting to note that two formula units (Ga2Se2) satisfy the 
18 electron rule (6 from 2 Ga and 12 from 2 Se) which is conducive to gap formation as 
in half-Heusler systems js]. The unit cell of /3 — GaSe however consists of 4 formula units 
containing 2 Ga dimers. 

Ab initio electronic structure calculations were carried out using the projector-augmented- 
wave (PAW) method, within density functional theory (DFT) as implemented in the 
Vienna Ab-initio Simulation Package (VASP) j?, 8]. The exchange-correlation potential was 
approximated by the Ceperley-Adler local density approximation (LDA) 9|]. This exchange- 
correlation potential was chosen over the gradient corrected version (GGA) because it is 
known that GGA underestimates the binding energies, which results in an overestimation 
of the lattice parameters lOj]. Since the interaction between the blocks of GaSe is week 
VdW type, the "GGA effect" becomes much more significant in the direction perpendicular 
to the atomic layers, resulting in a theoretical structure that is overly elongated along the 
crystallographic c-axis [4]. Such an elongated structure would not be appropriate for a 
theoretical investigation of the elastic properties of GaSe. In all calculations the outer s, p, 
d orbitals of the Ga and In atoms as well as the s and p orbitals of the Se were included 
in the valence states, while the rest were treated as core states. The cut-off energy for the 
plane wave basis was set to 300 eV and the convergence of self-consistent cycles was assumed 
when the energy difference between them was less than 10~^ eV. 

Various defect properties were investigated using a 3x3x1 supercell model containing 
72 atoms (36 Ga and 36 Se) (see Fig. [T](a)). Two types of In defects were studied, (i) 
In substituting for Ga (Inca) which disrupts the Ga-Ga dimer bond, and (ii) In at the 
interstitial site lui. Different charge states of the defects were investigated as a function of 
the Fermi energy. To identify the preferred location of the impurities, formation energies 
of the defects were calculated. All the atomic positions in the pure and defect containing 
supercells were fully relaxed until the quantum mechanical forces were less than 0.02 eV/A. 
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FIG. 1: (a) The crystal structure of layered GaSe. The interstitial site is located between the GaSe 
layers, equally far from the 6 NN Se atoms; (b) Formation energies associated with luQa and In;. 
When Ep is close to VBM (p-type GaSe), the In interstitial impurity becomes stable. 



The formation energy of a defect D in a charge state q denoted as D'^ is given by ll| : 

AHf {D^) = AE {D^) + ^ (rii/i,) + qEp 



(1) 



AE {D") = E {D'^) - E (GaSe) + {n,E{i)) + g^vBM 



(2) 



In eqs. ([T]) and ([2]), E {D'^) and E (GaSe) are the total energies of the defect-containing and 
the defect-free supercells and i?(i)'s {i = Ga, Se, defect) are the energies of the constituents 
in their standard solid state. The atomic chemical potentials /ij's are referenced to E[i) 
and nj's are the number of atoms removed from {rii > 0) or added to (n^ < 0) the system. 
The electron chemical potential (Fermi energy) Ep is referenced to the energy (ii^vBM) of 
the valence band maximum (VBM). In the present calculations (-Evbm) was determined as 
the average of the one-electron energy level of the VBM over the A;-points where the total 
energy was calculated. As pointed out by Zhang HI], this approach has the advantage 
that the band edges calculated this way are consistent with the defect transition levels and 
gives a better single-particle position for the shallow defects. Furthermore, the band gap 
of GaSe calculated with this "average band-edge" approach (E'^^'^^'^ = 1.68 eV^ is closer 
to the experimental value (E|^p- = 2.13eVj than the direct gap located at the F — point 
(Er^p = 0.85eV). 

To check the accuracy of our calculations vis-a-vis experiment we also calculated the 
charge transition level energies which correspond to the values oi Ep where the formation 
energies of a defect in two different charge states (g and q') are equal i.e. 
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e{q/q') = [AE{D'^') - AE{D'')]/{q - q') (3) 

Equation ([1]) shows that the formation energies of the defects depend on the atomic 
chemical potential (/i) of the constituents as well as on the charge state (g) of the impurity. 
The values of the /ij's are constrained by several physical conditions. For example: (a) to 
avoid precipitations, /ij's must be negative and (b) to maintain a stable host compound, the 
chemical potentials must satisfy /iGa + AtSe = AH (GaSe), where AH (GaSe) is the formation 
enthalpy of GaSe. The theoretical value of Ai7 (GaSe) = —1.12 eV. To avoid secondary 
phase formation between the host elements and impurities, one must impose additional 



constraints on the chemical potentials of the defects [ll|. However, the effect of these 
conditions would be a constant shift in the formation energies and since we are interested in 
the relative formation energies associated with different locations and charge states of the 
same defect we set ^def = /^in = 0. 

In Fig. [11(b) we give the calculated formation energies as a function of Ep for the 
two types of defects each in three different charge states. Also we give both theoretical 
{Ef;^' = 1.68 eV) and experimental (e^^p = 2 .13 eV^ band gaps. We find that he forma- 
tion energy of Itiq^ is 0.28 eV and one of the charge transition level associated with this 
defect, e(+l/0) = 0.28 eV above the VBM. This value is in fairly good agreement with the 
acceptor level at 0.21 eV, measured by Cui et al. using deep level transient spectroscopy 
(DLTS) |l^. This gives us confidence in our total energy calculations using DFT and the 
supercell model to understand the defect physics. 

The defect states with lowest formation energies are: Inf+ for Evbm(= OeV) < Ep < 
0.135 eV, In^t for 0.135 eV < Ep < 0.28 eV and In^^ for 0.28 eV < Ep. In^^ is the most 
stable defect for a wide range of Ep. However when the Fermi energy is tuned towards 
the VBM energy, Inf^ defect becomes most stable. We will discuss the underlying physics 
of this change by examining the single particle density of states (see below). As regards 
the effect of In^^ defect on the electronic structure, we find that the band structures near 
VBM and CBM are affected very little. One therefore does not expect much change in the 
transport properties in In doped GaSe if the impurity goes to a Ga site in the neutral charge 
state. One can understand this lack of significant change by looking at the In 5s — Ga4s 
dimer antibonding state (which hybridizes with the Se p-bands to give rise to states in the 
neighborhood of the band gap) and observe that it is not significantly different from the 




FIG. 2: (a) The total DOS of GaSe with Irij and the projected DOS of the In s-orbital, showing the 
positions of the HDDS (-5.5 eV) and DDS (just above Ep) introduced by the charged Ini defect; (b) 
The energy barrier which must be overcome in order to cleave the GaSe crystal increases dramat- 
ically when In occupies the interstitial site compared to the case when In occupies substitutional 
site. For comparison the case of pure GaSe is also shown. 

Ga4s — Ga4s dimer antibonding state. 

To understand why Inf^ has the lowest formation energy whenEp < 0.135 eV, we look at 
the electronic structure, the single particle density of states (DOS) and the nature of defect 
state introduced by Ini. Fig. [2](a) gives the total and partial (associated with In s) DOS for 
this case. We see that Ini introduces a hyper deep defect state (HDDS) near the bottom of 
the Se p bands (at ~ —5.5 eV). It is a bonding state formed out of In 5s and neighboring Se 
p states. The corresponding antibonding state splits off from the Se p valence band states 
and is denoted as the deep defect state (DDS). This picture is very close to what happens 
when In is a substitutional defect in PbTe 13|. The strong mixing between In 5s and the 
neighboring Se p— states leads to the removal of one state (per spin) from the Se p band 
which becomes the DDS. In terms of electron counting, two of the three electrons from In 
occupy the HDDS and the three electrons (two from the electrons occupying the valence 
band in pure GaSe and one from In) fill the DDS and partially occupy the conduction band. 
Thus lui acts like a donor. Since the three electrons occupy states with energies larger than 
EvBM; clearly the formation energy of In^ in charge state g = 0, 1, and 2 are higher than 
Inca for which neither the band structure nor the electron count change. By remo ving three 



electrons from Ini to obtain g = 3 charge state we can lower its formation energy 14\. 

Now that we understand the microscopic nature of the In defect states corresponding 
to two different positions we have to see whether the changes in the local bonding lead 
to appreciable modifications of the elastic stiffness of the GaSe matrix as seen experimen- 
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TABLE I: Elastic constants of Ga1_i.In2.Se. All values are given in GPa. 
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tally . To simplify our analysis of elastic stiffness we have used two approaches, one by 
calculating the elastic constants and the other by exploring the energy barrier involved in 
relative shearing of two blocks in a super cell. For Ini we use the second approach whereas 
for Inca (Gai_2;In2;Se) we use both the methods. Gai-^^InajSe is characterized by 5 elastic 
constants: Cu, C12, C13, C33, and C44. We have determined these quantities from total 



energy calculations for five different strain configurations [15|. For each x we have calcu- 
lated the theoretical crystal structures by minimizing the total energies with respect to the 
lattice constants: first with respect to the volume of the unit cell keeping the c/a ratio fixed 
and then with respect to c/a keeping the previously obtained equilibrium volume constant. 
The elastic constants for x = 0, 0.25, 1 were obtained using small unit cells (8 atoms/cell) 
and the Brillouin zone (BZ) was sampled by a F — centerd 12x12x3 fc— mesh. In the case of 
X = 0.0625 the calculations were performed on 2x2x1 supercell with the BZ sampled by a 
6x6x3 grid of fc— points. 

We observe a monotonic increase in the lattice constants as the concentration of the 
impurities (x-value) increases following the usual Vegard's law. For the end compounds GaSe 
and InSe the theoretical lattice constants are less than 3% smaller compared to experiment, 
while the c/a ratios are within 0.7% of the experimental values. The underestimation of the 
lattice parameters is due to the well-known overbinding effect of LDA. The calculated elastic 
constants are given in Table [H along with previously calculated theoretical results [l^ and 



experimental data 



171]. There is a good overall agreement with the available experimental 
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and theoretical values. The changes in the elastic constants of (Gai-^Jn^^Se) depend almost 
linearly on x: on one hand, we observe a monotonic decrease in Cn and C12, but on the other 
hand, the elastic constants C13, C33, and C44 show small enhancement with increasing x. 
This indicates that when In impurity occupies Ga sites, the crystal becomes softer in the a— 
and 6— directions (parallel to the atomic layers) and stiffer along the c— axis (perpendicular 
to the atomic layers). 

The effect of substitutional In on the elastic constants of GaSe can be understood if 
we examine the connection between the structural and electronic properties of the end 
compounds: GaSe and InSe. In the case of InSe the interlayer separation is 2.96 A, which 
is smaller compared to the interlayer distance in GaSe ( 3.15 A), indicating a stronger 
InSe interlayer interaction. Since the elasticity of the in the c— direction is determined 
predominantly by the interaction between the atomic layers, the crystal becomes stiffer with 
the increase of In concentration in GaSe. The softening of (Gai_a;Ina.Se) in the a— and 
6— directions with the increase in the composition x, can be easily understood, because the 
intralayer distances are longer and therefore the intralayer covalent bonds are weaker in InSe 
than in GaSe |18l . 1191 ] . Given that in the a- and b-directions there are no "interlayer regions" 
which could counteract the weakening of the atomic bonds, the crystal becomes softer as 
the In concentration increases. Although the substitutional In impurity seems to increase 
the elastic stiffness of GaSe along c— axis the, the effect is rather small (e.g. C33 increases 
by 7.8% from GaSe to InSe). So we have to look at a different mechanism for In induced 
increased inter-layer rigidity as seen experimentally. 

As we see the elastic constants do not change appreciably in In doped systems when In 
goes as a substitutional impurity. This is also seen in the calculations of energy barrier 
associated with relative shearing of two atomic blocks (each block being made up from 
4-atomic planes) in a unit cell. Fig. EJ^b) compares the energy barriers involved in this 
relative shearing. The energy barriers for both pure GaSe and Gai_a;Ina;.Se are very small 
and comparable. Thus substitutional In does not enhance the shear rigidity of GaSe. In 
the same figure we show the energy barrier associated with similar shearing in the presence 
of an interstitial charged In defect (Inf^). The energy barrier and the initial slope increase 
dramatically (by factors of ~ 10 and ~ 7 respectively) in the presence of Inf^. Clearly GaSe 
is very soft and inserting interlayer charged In defects can make the crystal rigid against 
shear distortion. We note that for certain shear configurations the interstitial In could not be 
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accommodated by atomic relaxation. Thus, for consistency, we present the results obtained 
for rigid shear. 

In summary, using ab initio electronic structure studies we have shown that In defects 
can go into GaSe in three ways, (i) mostly as a neutral substitutional defect, (ii) as a 
singly charged substitutional defect and (iii) as a triply charged interstitial defect Inf"''. The 
charged defects can be stabilized by p-doping GaSe. As a substitutional defect In does not 
impact appreciably the physical properties of the crystal. In contrast the interstitial In if in 
a proper charge state can increase the shear stiffness of GaSe [sl significantly. We suggest 
that the observed dramatic improvement when In is added to p-type GaSe is due to charged 
interstitial defects. One of the important predictions of our calculation is the presence of a 
deep defect state just above the valence band maximum when In is present as an interstitial 
defect. Deep level spectroscopy should be able to see this defect. 
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